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The rate for the lepton-number-violating decay K^ — ► tt^ // ± /U ± is calcu- 

q^ \ lated in a model which incorporates doubly-charged Higgs bosons. We find 
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that for reasonable values of the parameters the decay branching ratio may 
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D | be as large as 2 x 10 . Although this is a discouragingly small number, it 
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is of the same order of magnitude as the rate mediated by massive Major ana 

neutrinos. 
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Qh| In the past couple of decades a large amount of effort has been dedicated to the study of 

double beta decay in nuclei. The motivation has been primarily two-fold: The experimental 
observation and analysis of the two-neutrino process (A, Z) — » (A, Z + 2) + 2e + 2v has 

•i— I , 

shed light on the structure of the nuclear matrix elements involved, and the search for the 
neutrinoless process (A, Z) -^ (A, Z + 2) + 2e has provided a powerful tool for the study of 
the Majorana properties of neutrinos, with considerable speculation about what might be 
inferred from the actual observation of this lepton-violating decay. 

The studies of neutrinoless double beta decay have mostly relied on the assumption that 
the process would be mediated by massive Majorana neutrinos, with the modern picture 
being described by leptons which appear in left- and right-handed doublets, the large mass 
of the right-handed neutrino being responsible for the small left-handed neutrino mass via the 
seesaw mechanism. In this scheme neutrino double beta decay would proceed as a second- 
order process, described by the same Hamiltonian that is responsible for the standard weak 
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processes. 

However, the possibility that double-beta decay might be mediated by a different lepton- 
violating interaction was also considered early on |I[, with the suggestion that a A 3 / 2 ,3/2 
resonance inside the nucleus might induce a lepton-violating process. Later, a mechanism 
for |AL| = 2 was proposed where the electrons might couple to quarks via a doubly-charged 
Higgs @ , but subsequent arguments showed that in this case the effective coupling is negli- 
gibly small |§. In fact, the conclusion for nuclear double beta decay was that there are no 
important \AL = 2| contributions not proportional to m v in a SU(2) x U(l) model with a 
Higgs triplet. 

The Majorana nature of neutrinos can in principle also be tested in the double beta decay 
of strange mesons, although one would expect very small branching ratios. Experimental 
searches for nuclear double beta decay put stringent limits on K^ — > ir^e e . Calcula- 
tions for K^ —y 7t T ix^fx^ have been performed in a model where lepton-number violation is 
transmitted through the propagation of a heavy Majorana lepton [|]], with resulting rates of 
the order 10~ 16 sec -1 . A computation has halso been done using a relativistic quark model 
for mesons |5[], which made it possible to consider small neutrino masses. In that case, the 
resulting rates are several orders of magnitude smaller. 

Here we wish to look at the possible contribution to K — > -k^^^ that doubly charged 
Higgs bosons might make, and see how they compare to the prvious calculations just men- 
tioned. These doubly-charged Higgs bosons are present in exotic Higgs representations and 
can have lepton-number-violating couplings to any pair of leptons. They appear in many 
models, and to provide some focus we shall concentrate on the phenomenology of a model 
with one \Y\ = and one \Y\ = 2 triplet 0. The features of this model that are relevant 
for our analysis are the following: 

The lepton-number-violating coupling to left-handed leptons is specified by the following 
Lagrangian: 

L = ih ij ip? L CT 2 A4> jL + /i.e. (1) 



where ipn is the two-component leptonic doublet and the A is the 2x2 representation of the 
Y = 2 complex triplet. These couplings lead directly to Majorana masses for the neutrinos, 
so this model does not rely on right-handed lepton partners in order to generate neutrino 
masses. The Lagrangian in eq. (1) implies a Majorana mass term for the neutrino given by 

m „ = a*™. ( 2) 

where the subscript ij refers to the lepton families, v is related to the neutral Higgs vacuum 
expectation values, and sh is the sine of a doublet-triplet mixing angle. 

For the case of the electron neutrino, the experimental observation of nuclear double 
beta decay puts severe limits on the value of m^ and as a result on h ee . One obtains f| 

s H h ee < 6 x 1(T 12 . (3) 

Clearly for the electron-electron case we would not expect the model to predict any significant 
phenomenological consequences. However, in general for some models [0] there are certain 
niceties that follow from the assumption that the expectation value of the neutral Higgs 
vanishes. In that case the limit on m Ue is not relevant, and from Bhabba scattering one 
may derive the limit H 0] @ 

h\ e < 10~ 5 GeV~ 2 x m 2 A++ . (4) 

Constraints on h^ can be obtained from the experimental limit on muonium to antimuonium 
conversion, with the result 

heeh^ < 6 x 10- 5 GeV~ 2 x m\ ++ . (5) 

The physical states can be classified according to their transformation properties, and 
one finds a fiveplet H 5 (which includes the doubly-charged components), a threeplet H 3 , 
and two singlets. There are a variety of couplings of the Higgs bosons to fermions and to 
gauge vector bosons. The ones in which we are potentially interested are 



H£ + W W ~ iV2gm w s H g^ u 



H+ + H^W- ~ -iV2c H e(p ++ - p_) M (6) 

# 3 Vg ~ - ^ H [m q >{l + 75) - m,(l - 75)]. 

Finally, we remark that in the limit s# — > the standard model is regained, with the 
singlet H® playing the role of the standard Higgs boson with the standard couplings to 
fermions and vector mesons. 

II. K+ -> 7T-^V + 

The quark-level diagrams that contribute to the decay amplitude are shown in Fig. 1. 
For this particular model only the phenomenology of left-handed fields was developed M, 
but in any case limits on right-handed currents would make them irrelevant here. We adopt 
the usual form for the hadronic charged weak current 

J a = cos6 c [u ul a (l - 75K + «c7 Q (l - 7 5 K] 

+ sinfl c [M u7 Q (l - 7 5 K + «c7 a (l - 7sK] (7) 

The amplitude for Fig. 1(a) is then given by 



T la = — ■= < 7r| J^IO > — K-iV2gm w s H — T 
2V2 rriw m z 
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xh uu < lentous > —5 ^= < 01 J^\k > (8) 

m m^2V2 ' ' K J 

where < leptons >= u(h)(l + 7 5 )w c (/2) — (^1 <-► h)- As usual, the current matrix elements 
can be expressed as < 0| J^\k >= ifkPk and < n \J^\ >= if n p n , where f n and fk are the 
meson decay constants. The amplitude then reduces to 

T la = O.QG^mkfnfkE,, < leptons > H 2 — ■ (9) 

The amplitude for Fig. 1 (b) cannot be written down in a similar simple form. In order 
to estimate its contribution we begin by writing down the amplitude at the quark level: 



A = — — cos9 c u d ^ a (l - j 5 )u u ^-iV2gm w SH— — (10) 

xh UL , < lentous > — ^ = sin6U; s 7 Q (l — j 5 )v u . 

w m 2 w 2V2 



We project out the " pseudoscalar" content of A with the operation |TT)| 



Ti 6 = C n C k ^2 u u j 5 v s Av u j 5 u d (11) 

where the sum is over the quark spins and C n and Ck are normalization constants. The 
amplitude is then evaluated using a simple static constituent quark model. The decay 
constants are determined by applying the procedure to the decays K — > Iv and it — > Iv. 
The result yields Ck = ifkmk/8m s m u and C n = i/ 7r m, r /8m u m.d. The method can be tested 
by applying it to the amplitude of Fig. 1(a) and calculating the decay rate for that amplitude 
alone. In that case the result is the same as that found from the earlier calculation for the 
amplitude given in Eq. (0). 

Applying this procedure to Eq. ( |10"D yields 

T lb = -1.2 x 10 3 G 3 F /2 aC fc mV fc K < leptons > %-^-. (12) 

m 2 A++ 

The remaining diagrams in Fig. 1 can be evaluated in similar fashion. For example, T\ e 

is given by 



Ti e = i t=- sin Q w sin 9 C cos 9 C — ^ — 5 5 — < leptons > 

2V2 m% / m\+m\ ++ 

x(m s - m u )v«{l + 7 5 )^n d 7 a (l - 7 5 K(p++ +P-T- (13) 



where p ++ is the momentum of the incoming A ++ and p~ the momentum of the outgoing 
A_. 

It is clear that with reasonable values for the Higgs masses the contributions from these 
graphs will be considerably smaller than that from the first two. Therefore we calculate 
the rate using only the diagrams in Fig l(a,b). After squaring the amplitude and doing the 
spin sum we integrate over phase space as follows: 



Here s = (pi + P2) 2 , Pi and P2 are the muon momenta, and as usual X(x, y, z) = x 2 + y 2 + 
z 2 — Ixy — 2xz — 2yz. Evaluation of the integral yields the following result for the rate: 

( sijhfji, \ „ _ .a 1 



i? = 6xl0- 2 (Z^m) GeV 4 — . (15) 

\m z A++ J sec 



III. RESULTS AND DISCUSSION 

It is possible to develop models that incorporate Higgs triplets where the 
A ++ W~W~ coupling is allowed only with an enabling non-zero vacuum expectation value for 
the neutral member [fj]]. In that case one may obtain a limit for h ee through the Majorana 
mass term in neutrinoless nuclear double beta decay, similar to that shown in Eq. (J|). The 
assumption that the vev is zero avoids problems with maintaining p = 1 , but in that case 
graphs like Figs. l(a,b) would not be allowed. However, couplings of the triplet to vector 
bosons of the type shown in Eq. @ are a natural consequence of the model being discussed 
here. 

At the same time, the coupling of Higgs to fermions is determined by the overlap of the 
mass-eigenstate Higgs fields with the doublet, and the H 5 does not have such overlap. The 
H 3 does and that makes graphs in Figs, l(c-h) possible, but the contributions to the rate 
are smaller by at least a factor ~ (mk/m 1 \+) A . 

It is clear then that the rate will be dominated by Figs. l(a,b), with the result shown in 
Eq. fllBD and with a branching ratio given by 

BR = 0.9 x 1(T 9 ( S JlhBL\ G eV 4 . (16) 

\ m A++ / 



The existing experimental limit [11] for this branching ratio is 1.5 x 10 , corresponding to 



SHh m < 1.5 x 10 2 (m A++ /GeV) 2 . (17) 

This result is not invalidated by the previous limits in Eqs. (0fj,^), but it is not exactly a 
giant step forward either. Taking the values mA++ ~ 100 GeV |7| and ssh^ ~ 5x 10~ 2 mA++ 
(near to and not excluded by the limits of Eqs. (§,|5|)) yields a branching ratio of 2 x 10~ 16 , 
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which is discouragingly small but of the same order of magnitude as the rate for kaon double 
beta decay induced by light or heavy Majorana neutrinos and reported in earlier publications 

!§• 

In conclusion, we have calculated the rate for K ± —>■ -k^^h^ that can be expected for a 
decay mediated by doubly-charged Higgs. The branching ratio for reasonable values of the 
model's parameters is extremely small compared to the experimental limit, but of the same 
order of magnitude as that induced by virtual Majorana neutrinos. Since the question of 
lepton-number violation is still outstanding, these processes warrant further theoretical and 
experimental study 
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Figure Captions 
FIG. 1. Diagrams contributing to K + — > n~fx + fx + . 
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